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Characterising Body Tissue 



Field of the Invention 



The present invention relates to methods for the characterisation of body tissue. More 
specifically, the invention is concerned with the characterisation of body tissue as normal (e.g. 
healthy) or abnormal (e.g. pathological). The invention has particular, although not 
necessarily exclusive, applicability to the diagnosis and management of cancer, including 
breast cancer. 



Background 

In order to manage suspected or overt breast cancer, tissue is removed from the patient in 
the form of a biopsy specimen and subjected to expert analysis by a histopathologist. This 
information leads to the disease management program for that patient. The analysis requires 
careful preparation of tissue samples that are then analysed by microscopy for prognostic 
parameters such as tumour size, type and grade. An important parameter in tissue 
classification is quantifying the constituent components present in the sample. Interpretation 
of the histology requires expertise that can only be learnt over many years based on a 
qualitative analysis of the tissue sample, which is a process prone to intra observer variability. 

Despite the relative value of histopathological analysis, there remains a degree of imprecision 
in predicting tumour behaviour in the individual case. Additional techniques have the potential 
to fine-tune tissue characterisation to a greater degree than that currently used and hence will 
improve the targeted management of patients. 

< 

In existing research in this field, x-ray fluorescence (XRF) techniques have been used to 
study trace element composition of breast tissue and have shown that breast cancer is 
accompanied by changes in trace elements and such measurements could contribute to 
tissue grading \ It has also been shown that x-ray diffraction effects can operate as an 
effective means of distinguishing certain types of tissue 2 ' 3 . Furthermore, it has been shown 
that such diffraction effects could be suitably analysed to demonstrate small differences in 
tissue components and that this analysis could lead to a quantitative characterisation of 
tissues 4 . 

There remains in particular a need for techniques that can be used to characterise tissue to a 
greater degree in vivo, in order that the need for often painful and distressing biopsies can be 
reduced. 
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Summary of the Invention 

It is one general, preferred aim of the present invention to develop quantitative analytical 
approaches that could add precision to the characterisation of tissues, in particular to 
distinguish between normal and diseased (e.g. pathological) tissue. 

A preferred aim is to add precision to the several subjective components of tissue analysis, 
most notably those variables 'scored* in breast tumour grading. 

In general terms, the invention provides methods for analysing and/or characterising body 
tissue in which results are obtained by considering a combination of two or more different 
types of measured tissue properties. 

In one aspect, the present invention provides a method analysing body tissue, the method 
comprising: 

obtaining data representing a first measured tissue property of a body tissue sample; 
obtaining data representing a second, different tissue property of the tissue sample; 

and 

using the data in combination to provide an analysis of the tissue sample. 

In another aspect, the invention provides a method for characterising body tissue, the method 
comprising: 

obtaining data representing a first measured tissue property of a tissue sample; 
obtaining data representing a second, different tissue property of the tissue sample; 

and 

using the data in combination to provide a characterisation of the tissue sample. 

The characterisation in the second aspect may consist of characterising the tissue sample as 
normal or abnormal. Alternatively, the characterisation may be performed accounting for 
many grades of abnormality, for example, on a scale with "Normal" at one end and 
"Abnormal" at the other, with numerous positions therebetween. As a further alternative, or 
perhaps additionally, the characterisation may take the form of tissue typing, wherein the 
characterisation includes an expression of a specific trait, such as the kind of tissue, or the 
stage of cancer and the like. 

In either aspect it is preferred that data representing a third measured tissue property is also 
used in combination with the other data in the analysis or characterisation of the tissue 
sample. 

It is particularly preferred that data representing four or more measured tissue properties is 
used in combination in the analysis or characterisation of the tissue sample. 
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Suitable techniques that can be used to obtain the tissue property data include x-ray 
fluorescence (XRF), energy or angular dispersive x-ray diffraction (EDXRD), Compton scatter 
densitometry, low angle x-ray scattering and the measurement of linear attenuation 
(transmission) coefficients. 

The tissues properties that are measured may include the composition of the tissue sample, 
for instance the presence, concentrations and/or proportions of specific elements or organic 
compounds. Indeed, a tissue sample may contain more than one type of tissue, for example 
- fatty or glandular for instance, and the measured property may include information relating 
to this. 

Preferably, in either of the aspects above, the data is used in combination to obtain the 
desired result by using the data as the input to a predefined calibration model that relates the 
combined data to one or more tissue characteristics (e.g. normal or abnormal). 

In a further aspect, the invention provides a method for creating a tool for the analysis and/or 
characterisation of body tissue, the method comprising creating a calibration model that 
relates data representing two or more (preferably three or four or more) measurable tissue 
properties to one or more tissue characteristics. 

The calibration model is preferably produced by using sets of the measured data from tissue 
samples for which the characteristic(s) (e.g. normal / abnormal) to be determined by the 
model are already known. These data sets can be used to 'train' the model in a known 
manner. 

Other multivariate analysis techniques may be employed. 

It is a further general aim of the present invention to provide a method of analysing and/or 
characterising based on a recognition that Compton scattering densitometry techniques can 
be used in the analysis of body tissue to very effectively discriminate healthy and abnormal or 
diseased tissue and to discriminate types of abnormal tissue. Moreover, Compton scattering 
has been recognised as having potential for application to In vivo tissue characterisation 
techniques. 

The invention of this further general aim provides a method for analysing and/or 
characterising body tissue, the method comprising: 

obtaining Compton scatter data measured from a body tissue sample on which a 
penetrating (e.g. X-ray) radiation beam is incident; and 

using the data to provide an analysis and/or characterisation of the tissue sample. 

Compton scatter results from an interaction that occurs between a photon and an electron. 
For this interaction the electron is assumed to be unbound and acting as a free particle. This 
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assumption can be made if the energy of the incident photon is much greater then the binding 
energy of the atom. Figure 8 illustrates the Compton interaction, where E 0 is the energy of the 
incident photon, Ei is the energy of the scattered photon, m 0 c 2 is the rest mass energy of the 
electron and 9 is the scattering angle of the photon and cp is the scattering angle of the 
electron. T is the kinetic energy imparted to the electron. 

The electron taking part in the interaction is assumed to be stationary, i.e. the initial energy 
(E e ) and momentum of the electron equals zero. During the interaction the photon imparts 
some of its energy to the electron. The amount of energy transferred determines the angle of 
the recoil of the electron and the angle of the resultant photon. 

The angle and energy of a Compton scattered particle can be accurately calculated using the 
principle of conservation of energy and momentum. From Figure 8 it can be seen that the 

incident photon has energy E 0 =hv and the scattered photon has energy -hv\ 

Resolving the energy and momentum into parallel and perpendicular components gives the 
important Compton scatter equation 
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hence a measure of Compton scatter can be made by detecting the appropriate energy 
photons at a given angle. 

In some instr?; ces it may be sufficient for the Compton scatter data to be as simple as a count 
of photons detected at a selected angle/energy in a given time period. In other instances, it 
may be desirable to obtain an absolute measure of electron density (or some other derived 
measurement). Particularly in the latter case, the Compton scatter data is preferably 
corrected for attenuation in the tissue sample. 

One way to compensate for attenuation effects is to use two radiation sources and two 
detectors. This is an approach commonly used in bone densitometry, but is less preferable 
when examining tissue samples, particularly in vivo, because it results in a greater dose of 
radiation. 

A preferred method to correct for attenuation effects is to obtain data representing a measure 
of the directly transmitted x-ray radiation for each Compton scatter measurement. This data 
can then be used to correct the Compton scatter data for attenuation in the tissue sample. 

Especially at low angles (less than 90°), it is also important to be able to distinguish Compton 
scatter measurement from the coherent scatter peak. So, where the transmitted radiation is 
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to be used to correct for attenuation it is preferable that the energy of the scattered photons 
detected is as close as possible to that of the transmitted radiation. This ensures that the 
attenuation coefficients are not too different for the two measurements. The energy of the 
incident penetrating radiation beam and the angle selected for Compton scatter measurement 
are chosen such that the Compton and coherent scatter peaks can be resolved, whilst 
minimising the separation (i.e. energy) of these peaks. This substantially eliminates self- 
attenuatlon effects as it allows one to assume that the attenuation coefficients in the sample 
affecting both peaks are substantially the same. 

Preferably the data is used as the input to a predefined calibration model that relates the 
Compton scatter data to one or more tissue characteristics (e.g. normal or abnormal or a 
scale of abnormality having "normal" at one end of the scale and "abnormal" at the other). It 
is particularly preferred that the Compton scatter data is used as an input to a multivariate 
model. 

Although the specific embodiments describe the various aspects of the invention in relation to 
breast cancer, it is to be understood that the invention, generally, has a much wider 
applicability. Indeed, along with analysing and characterising breast tissue for cancer other 
assessments, such as general nodal assessment, liver, pancreas, prostate, colorectal 
assessments are contemplated, also urological and gynaecological assessments are also 
envisaged. 

Brief Description of the Drawings 

Embodiments of the invention are described below by way of example with reference to the 
accompanying drawings, in which: 

Figure 1 is a schematic diagram of EDXRD experimental apparatus employed in the 
exemplary methods described below according to embodiments of the invention; 

Figure 2 is a series of graphs showing average XRF responses; 

Figure 3 shows EDXRD scatter profiles for normal and diseased tissue; 

Figure 4 shows PLS model predictions for the normal test samples; 

Figure 5 shows PLS model predictions for the diseased test samples; 

Figure 6 shows predictions of tissue type for the normal test samples; and 

Figure 7 shows predictions of tissue type for the disease test samples; 

Figure 8 illustrates the energetics of Compton scattering; 
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Figure 9 shows schematically the experimental set-up used for Compton scatter 
measurements in an example of an embodiment of the invention; 

Figure 10 illustrates the sample holder used in the Compton scatter measurement of the 
example; 

Figure 1 1 shows the peak measured with the Ortec GLP-25300 HPGe detector, used in the 
experiment, using an Am-241 source; 

Figure 12 is a schematic of the electronics used for electron density measurements; 

Figure 13 shows an observed scatter spectrum obtained for one sample during the 
experiment; 

Figure 14 shows the apparatus of Figure 2 set-up to take transmission measurements; 

Figure 1 5 shows a calibration graph for the electron density measurements; 

Figure 16 is a graph of differential scatter coefficient versus theoretical electron density; 

Figure 17 shows the results from the Compton scatter measurements taken from all samples 
during the experiment; 

Figure 18 is a graph of tabulated tissue values and experimental data. 

Figure 19 illustrates the cylindrical geometry used as the sample holder for the measurement 
of the electron density; 

Figure 20 shows a scatter spectrum from a malignant breast tissue sample; 

Figure 21 shows a calibration graph of the calculated linear scatter coefficients against the 
counts measured in the Compton scatter peak for the calibration solutions; 

Figure 22 shows a graph of the differential scatter coefficient from experimental data against 
the calculated electron density for the calibration solutions; 

Figure 23 shows a box plot of the electron density results obtained from the tissue samples; 
and 

Figure 24 shows a graph of the electron density values for each tissue type. 
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Description of Embodiments 

The following describes a study that has been conducted in vitro using breast tissue samples. 
The principles exemplified below are, however, more widely applicable, for example to other 
tissues and in relation to data obtained from in vivo as well as in vitro measurements. The 
principles can usefully be employed in in vivo imaging applications. 

The exemplary embodiments described below employ x-ray fluorescence (XRF) and energy 
dispersive x-ray diffraction (EDXRD) techniques to reveal the tissue characteristics. The 
invention is not, however, limited to these two techniques and other techniques may be used 
in addition or as alternatives to XRD and EDXRD. Other techniques that might be used 
include Compton scatter densitometry, low angle x-ray scattering and linear attenuation 
(transmission) coefficients. 

As discussed in more detail below, during the study the concentrations of K, Fe, Cu and Zn 
were measured in 77 breast tissue samples (38 classified as normal and 39 classified as 
diseased) using X-Ray Fluorescence (XRF) techniques (in other embodiments, 
concentrations of other elements or organic compounds might be measured). The coherent 
scattering profiles were also measured using Energy Dispersive X-Ray Diffraction (EDXRD), 
from which the proportions of adipose and fibrous tissue in the samples were estimated. 

The data from 30 normal samples and 30 diseased samples were used as a training set to 
construct two calibration models, one using a Partial Least Squares (PLS) regression and one 
using a Principal Component Analysis (PCA) for a Soft Independent Modelling of Class 
Analogy (SIMCA) technique. The data from the remaining samples, 8 normal and 9 diseased, 
were presented to each model and predictions were made of the tissue characteristics. 

Three data groups were tested: XRF. EDXRD and a combination of both. The XRF data 
alone proved to be most unreliable indicator of disease state with both types of analysis, il.e 
EDXRD data was an improvement, however with both methods of modelling, the ability to 
predict the tissue type most accurately was by using a combination of the data. 

1. Breast Tissue Samples 

The tissue samples measured were obtained from mastectomies, lumpectomies and breast 
reduction surgery. In regard to the latter, a number of healthy breast tissue samples were 
obtained. The tissue obtained from mastectomies or lumpectomies, was generally taken from 
the site of a lesion, classified as invasive ductal carcinoma, and in some cases normal tissue 
was taken from areas distant to the tumours, in line with the available samples, 
investigations were made for 38 samples classified as normal and 39 samples classified as 
diseased. The weight of each of the specimens was of the order of 1 g. Most specimens 
were of thickness in the range of 2-3 mm. Following excision the samples were kept frozen at 
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- 85 C, no processing or sample preparation taking place between excision and 
measurement. For both the XRF and !h3 EDXRD measurements the samples were allowed 
to thaw before being measured in room temperature. 

2. Experimental Procedure 
2.1 XRF 

The XRF studies were carried out making use of the European Synchrotron Radiation Facility 
(ESRF), working on the Bending Magnet beamline BM28 5 . Using a simple arrangement of 
incident synchrotron radiation, tuned to a photon-energy just above the K-absorption edge of 
interest, particularly low elemental detection limits are achievable (< 1ppm). The high 
intensities of XRF available allow for short measurement times, providing for a high sample 
throughput. For the synchrotron photon beam, the plane of polarisation is the same as that of 
the electron orbit. Thus for a 90° geometry between photons directed on to the sample and 
the normal to the detector (Si(Li), Gresham Scientific Instruments, Sirius model), the strong 
linear polarisation of the photon beam provides significant suppression of the scattered 
photon intensity (fluorescence being unaffected). Given that the detector has lateral extent, 
the remaining sample-dependent scattered photon (coherent and incoherent) intensity 
reaching the detector is therefore governed by the solid angle formed between the sample 
and detector crystal. In addition to providing improvement in the signal to background ratio, 
control of the scattered radiation intensity allows use of the scattered peak area as a 
normalisation factor. As tissue is a low Z material, the fact that the detection system cannot 
resolve the Compton component will not affect the results. 

Each element of interest (K, Fe, Cu and Zn) was identified by the photopeak associated with 
its K* fluorescence photon emission. In order to seek maximum production of the K« photons 
of interest, the samples were irradiated by photons of energy 500 eV above the particular K 
absorption edge, being an arrangement which also allows for the resolution of the scattered 
incident peak and the fluorescence response. The exception to this method was for K, where 
the data were collected using the same incident photons as that for Fe. In order to quantify 
the sample concentrations of the elements of interest, calibration curves were constructed for 
each element. 

The calibration standards were aqueous solutions of the elements, the water matrix of the 
calibration models matching the "wet" nature of the tissue specimens. The following ranges 
of concentrations were used for the calibrations, as indicative of those expected to be found in 
tissue: 

K: 100, 300 and 1000-4000 ppm in increments of 1000 ppm 
Fe: 3-30 ppm in increments of 3 ppm 

■ 
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Cu: 1-10 ppm in increments of 1 ppm 



Zn: 2-25 ppm in increments of 2 ppm 



The calibration solutions were measured in petri dishes that were sealed with laboratory 
sealing film (LabSeal, Merck). The tissue specimens were placed on such petri dishes 
previously filled with purified water and sealed. The specimens were then covered with the 
same sealing film. The beam size on the specimens was 3 mm x 0.5 mm. 

The spectra acquired from the standard solutions were analysed using the software PeakFit 
(PeakFit™ SPSS Inc. AISN Software Inc) developed for spectroscopy. The spectra were 
smoothed using a procedure based on deconvotution, leading to the removal of peak 
broadening effects caused by imperfect resolution of the measuring instruments. The spectra 
were subsequently fitted using a procedure based on the Levenburg-Marquardt non-linear 
minimisation algorithm. The fitting process took into account a linear baseline resulting in an 
estimation of the net total of counts integrated over the width of the photopeak. In order to 
normalise the fluorescence response, the scattered photopeak area was also calculated. The 
ratio of fluorescence to scattered photon peak area was then used to derive the relationship 
between element fluorescence and its concentration. 

The tissue samples were irradiated under the same conditions as those used for the standard 
solutions and spectra were collected for each of the elements of interest. The analysis of the 
spectra also followed the procedure described for the standard solutions. The least squares 
fit derived from the calibration data, which relates the ratio of fluorescence to scatter photon 
peak area and the element concentration was used to quantify the levels of each element in 
each of the samples. It is acknowledged that only a small area of the samnle is irradiated but 
measurements indicate that the inhomogeneity has been found not to significantly alter the 
profound differences between healthy and cancerous tissues. No correction has been made 
for matrix effects has been made in this study. However, the interest here is in the 
comparison between the levels of healthy and cancerous tissue and as any errors are 
systematic the comparison is not compromised. 

2.2 EDXRD 



The EDXRD scatter profile of each of the samples was measured using a technique that 
utilises the scatter of a polyenergetic photon beam at a fixed scatter angle. This technique 
has been used for a number of biomedical applications, notably that of estimating bone 
mineral density 6,7 and more recently for breast tissue analysis 2 . For an overview of 
applications of X-ray diffraction analysis in crystalline and amorphous body tissue see as ai 
instance Royle et al 8 . 



9 



WO 2005/055827 



PCT/GB2004/005185 



The experimental set up is shown schematically in figure 1 . A tungsten target x-ray tube 
operating at 70 kV and 1 5 mA was used, the intrinsic filtration being 1 mm of beryllium. The 
incident beam is tightly collimated via a slit cut in a dural slab forming a rectangular cross 
section of dimensions 1 mm x 2 mm. A similar collimation arrangement was set up at a 
scatter angle of 6° leading to a scattering volume in which the thickness of the sample was 
enclosed. The samples were mounted in standard 35 mm slide frames, sealed on either side 
with film 4 microns thick (Ultralene from Glen Spectra Reference Materials). The translator 
enabled the samples to be moved through the incident beam, in this instance a distance of 3 
mm, producing an approximate irradiated volume of 12 mm 3 . The scattered photons were 
detected using an HpGe detector (EG&G Ortec). The output pulses were analysed using a 
multi channel analyser (92X Spectrum Master, EG&G Ortec). The measurement time for 
each sample was 2400 seconds. 



3. 



Results 



3.1 



XRF 



Figure 2 shows the average of all XRF spectral response for the four elements of interest for 
normal and diseased tissue. Quantitative values of the elemental concentrations were 
obtained from the ratio of the XRF response peak to the scattered peak (not shown in figures) 



via the calibration line. 



Table 1 summarises these results for all tissue samples, showing the range of concentrations 
within a sample group, the mean concentrations and the standard deviation. For details 
about the measured levels of the elements and the comparison between healthy and 
cancerous tissue see work by Geraki etal\ 



Units in 


K 




Fe 


Cu 


Zn 


PPM 


Norma 
I 


Diseased 


Normal 


Disease 
d 


Normal 


Disease 
d 


Normal 


Disease 
d 


Range 


57- 
514 


114-1120 


< D.L- 
36.64 


< D.L.- 
48.98 


< D.L- 
1.53 


< D.L.- 
2.09 


< D.L- 
7.51 


< D.L- 
20.19 


Mean 


163 


512 


8.82 


17.62 


0.27 


0.88 


2.19 


6.96 


~SD 


"141 


336 


~8J4 


"t4!09 


0.22 


0.52 


1.37 


5.31 



Table 1 : Summary of the XRF results for the elements K, Fe. Cu and Zn. 



3.2 



EDXRD 



Figure 3 shows the averaged diffraction spectra for all normal and all diseased tissue 
samples. The difference in the composition of the two types of specimens is evident. 
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The characteristic peak from the adipose tissue can be seen at a momentum transfer value of 
1.1 nm and the characteristic peak from fibrous tissue is at approximately 1 .6 nm" 1 . These 
peaks were fitted using the same technique as for the XRF spectra. The evaluated 
photopeak areas reflecting the presence of adipose and fibrous tissue were then corrected for 
factors such as the shape of the x-ray tube spectrum and the difference in the attenuation and 
scatter properties of the two types of tissue. The corrected relative intensities of the two 
scatter peaks reflect the relative amounts of the two materials in each specimen. 

These results show a strong dependence of the relative amounts of the two materials on the 
type of specimen. The healthy specimens were predominately made of fat (76 ± 9 %) while 
the tumour specimens were mainly composed of fibrous tissue (85 ± 4 %). 

4. Modelling 

The above data were divided into two groups i.e. 30 normal and 30 diseased sample data 
were used as a training set to produce a calibration model. The data from the remaining 
samples from each group (8 normal and 9 diseased) were then used as input for the model 
and the tissue type was predicted. 

In the calibration process we use the empirical data (i.e. elemental composition and EDXRD 
data) and prior knowledge (i.e. quantities known to represent that data) to predict unknown 
quantitative information from future measurements. In this simple application we are using a 
multivariate approach i.e. we use many variable measurements X{x(1), x(2)....x(n)} to quantify 
the target variables Y{y(1), y(2)....y(m)}. In this case the X variables are the data measured 
above and the Y variables are the categories of tissue i.e. normal or diseased. 

In order to predict the later from the former we need to estimate how X relates to Y and an 
example of this is a regression model of the form y=Xb+c. The unknown parameters b and c 
are estimated from the calibration data which can then be used to predict future values of y 
from measured X. 

Another method is to use a classification procedure where models are created that represent 
a particular classification of variable. This is carried out using principal component analysis 
on the data belonging to a particular classification. Input data are then analysed and 
compared with the models and a fit to each classification is established. 

There are several statistical calibration methods available for creating a model 9 and this work 

uses a partial least squares (PLS) regression and principal component analysis (PCA), which 

is used for Soft Independent Modelling of Class Analogy (SIMCA) classification. The 

Unscrambler 10 is a software package that was used to perform multivariate analysis in this 
study. 
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4.1 Partial Least Squares (PLS) regression 

Three models were created using the PLS approach. The first used the XRF data alone, the 
second used the EDXRD data alone and the third used the combined data set. The 
appropriate data from the test samples (8 normal, 9 diseased) were then presented to the 
mode! and predictions of tissue type made. 

Figures 4 and 5 show the predictions for the normal samples for all three models and the 
predictions for the diseased samples respectively. 

4.2 Classification models 

For each of the data groups (XRF, EDXRD and combined data) principal component analysis 
models were constructed for the healthy and the cancerous samples. The appropriate data 
from the test samples were presented to the models and a score was obtained indicating how 
close each sample came to each of the models. 

Figure 6 and 7 show the predictions for each sample in the normal category and diseased 
category respectively. 



5. 



Discussion 



It can be seen that the XRF model predictions using the PLS approach are the most 
unreliable with the EDXRD and the combined data being similar. However the uncertainty in 
the predictions becomes significantly smaller when both data sets are combined. 

If an acceptable prediction parameter was chosen at 70% certainty and above, the number of 
reliable and false predictions can be found. Table 2 summarises the data for the PLS model 



Mean 
Prediction 
(%) 



Mean 
Uncertaint 

y(%) 



True +ve 
False -ve 



XRF 



Normal 



72.8 



11.2 



Diseased 



69.8 



32.9 



EDXRD 



Normal 



82.2 



20.1 




Diseased 



98.5 



21.6 



Combined 



Normal 



82.6 



12.2 



Diseased 



99.2 



14.5 



Table 2 : The mean predictions, uncertainty and number of true, false and undecided 
predictions for each data group and tissue type using PLS. 
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Similarly, when using the classification technique, the XRF is the most unreliable with an 
improvement being shown using the EOXRD data. It should be noted that for the normal 
samples predictions were higher in the normal classification whereas for the diseased 
samples the wrong classification was made in several instances. The use of the combined 
data shows a marked improvement in prediction particularly when examining the diseased 
samples. 

As above, choosing a 70% probability cut off limit, the number of accurate predictions can be 
found. Table 3 summarises the data. 




Table 3.: The mean predictions, true and false positives for each data group and tissue type 
using the classification technique. 

The relative inefficiency of the XRF data (compared to the combined XRF and EDXRD 
results) in successfully and accurately predicting the type of the test specimens is due to the 
wide spread of concentrations that characterise the groups of samples, evident by the large 
associated standard deviations (table 1 ). 

As illustrated by the study described above, embodiments of the present invention can 
provide improved characterisation of tissue types using a combination of data and an 
appropriate model. A classification technique has been shown to be particularly successful. 

Embodiments of the first general aim of the present invention have been described above by 
way of example. It will be appreciated that various modifications to that which has been 
specifically described can be made without departing from the invention. For instance, the 
study described above to exemplify the invention involved the use of only two types of tissue 
characterising properties. Other embodiments of the invention may use more than two types 
of tissue characterising properties or alternative characterising properties. Creating a model 
using samples that are characterised using a variety of useful parameters may develop useful 
histopathology tools. Provided the different data groups can represent all the parameters one 
wishes to quantify, the multivariate approach is a promising method for accurate 
characterisation of samples. 
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Although the specific embodiments described above relate primarily to breast cancer, it is to 
be understood that the invention, generally, has a much wider applicability. Indeed, along 
with analysing and characterising breast tissue for cancer other assessments, such as 
general nodal assessment, liver, pancreas, prostate, colorectal assessments are 
contemplated, also urological and gynaecological assessments are also envisaged. 

Compton Scattering 

The invention of the second general aim is exemplified below with reference to in vitro 
Compton scatter measurements from uniform samples of body tissue. The general technique 
is, however, equally applicable to the analysis of non-uniform tissue samples, including in vivo 
applications. 

The experiment was undertaken twice (A & B), each time comprising two sections; Compton 
scatter measurements were made on all the samples, followed by transmission 
measurements. This was done in preference to the two measurements being made 
consecutively for each sample. This method was adopted for two reasons; firstly to ensure 
consistency of set-up between samples through minimising the moving of equipment and 
secondly to save time. 

Theory 

The angle and energy of a Compton scattered particle can be accurately calculated using the 
principle of conservation of energy and momentum. If the incident photon has energy 
E, = h vand the scattered photon has energy E 2 =hv'. Resolving the energy and 

momentum into parallel and perpendicular components gives the important Compton 
equation 

*> =—7-^ (1) 



1 + 



m Q c 2 j 



(l-costf) 



where m 0 c 2 is the rest mass energy of the electron and 9 is the scattered angle. 
Consider the cylindrical geometry shown in figure 1 . A beam of photons is in the direction AB 
with energy E, and a detector is placed at a scatter angle 6 to the incident beam. The 
number of scattered photons, S, with energy E 2 reaching the detector is given by 

S oc (Vp e )exp(- £//, (x)dx j exp^- f Ml ( x )dx) ( 2 ) 

where V is the volume of scattering material. p e is the electron density of the material in the 
scattering volume, jj. is the attenuation coefficient of photons at the incident energy, and u? is 



14 



WO 2005/055827 PCT/GB2004/005185 



the attenuation coefficient of the Compton photons, scattered through angle 9 with reduced 
energy E 2 rif the incident energy E t and scatter angle 8 are carefully chosen it can be 
assumed that E^Ea and therefore that y,- n 2 . Usin 9 these assumptions it follows that 

and equation (2) becomes 

S «c(Fp,)exp(-f A (x>fa) (4) 
From the exponential law of attenuation it is found that 

— = exp (- J/i, (x)dx) (5) 



where l 0 is the incident photon intensity and I is the transmitted photon intensity. Therefore by 
obtaining a measure of i 0 , 1 and S, the electron density can be found from 



p e =k 



-1 



(6) 



where S is the scattered count intensity, T=l/I 0 and k is a constant which includes the volume 
term determineu using calibration materials with known electron densities. 

EXPERIMENT A 
Method 

Ideally a monoenergetic source should be used to ensure that the Compton scatter peak is 
easily detectable. The characteristic lines produced by the x-ray tube were used to generate a 
pseudo-monoenergetic source. Using this method the Compton and Coherent peaks can be 
easily resolved and windowed. The bremsstrahlung background can then be subtracted. 

The desired outcome from the experiment was to be able to resolve the Compton and 
coherent peaks, whilst keeping them as close in energy as possible. The detector 
characteristics dictate that the minimum resolvable energy is about 1 keV. 
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Using the Compton scatter equation set out further above the angle required to give an 
appropriate difference in energy between the incident and scattered peaks was calculated. To 
obtain Eq - E 1 s 1 keV, there was a choice between either using a higher energy and a small 
scattering angle (defined as 0 in Figure 9 below) or a lower energy with a larger scattering 
angle. Both of these options were explored. The conclusion was that using a higher energy 
and minimising the angle offered a number of advantages. Firstly the attenuation of the beam 
by tissue* will be lower with a higher energy. Secondly the scatter is at a maximum in the 
forward direction and at a minimum at 90°. Therefore the flux reaching the detector will be 
much higher with a smaller angle, reducing count times considerably. A smaller beam size 
can also be used, improving the accuracy of the measurement. Larger scattering angles can, 
however, be used if desired. 

The distance between the source, sample and detector were kept to a minimum to decrease 
the loss of flux due to inverse square law effects. The experimental set-up is shown in figure 
9. 

The incoming x-ray beam was collimated to a 0.5mm circle, both before and after the sample. 
This was the smallest beam size obtainable whilst maintaining a reasonable flux. The Ka lines 
from the tungsten target of the x-ray source (Ek^sg.SkeV and £^=57.97 keV) were used. 
At this energy a scattering angle of 30° would give a peak separation of 1 keV between the 
Compton and coherent peaks. The exact angle that was set-up in this example was 28.2°. 
The scattering volume comprises the tissue contained within the intersecting area of the 
incoming and scattered beam. For this beam collimation and a scattering angle of 28.2° the 
entire scattering volume was contained within the sample. This means that no air or plastic 
was contained within the scattering volume. The samples were measured for 20 minutes per 
position for 12 positions around the sample. 

Samples 

5 samples of each tissue type were chosen for examination. These were 5 fibroadenoma 
(benign), 5 invasive ductal carcinomas (malignant) and 5 pure adipose (normal) samples. The 
samples were placed into plastic pots of 6mm inside diameter and imm wall thickness 
(illustrated in Fig. 10). Although the walls of the container were relatively thick and would 
cause significant attenuation of the scattered beam, these containers were chosen because 
they offered a number of important advantages: 

1 ) The sides were completely rigid so the samples could be placed into the pots and 
lightly compressed with a stopper (see Figure 10) to remove any air gaps without the 
pot becoming distorted. This also minimises tissue movement throughout the 
experiment. 
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2) The pots were cheap so each sample could have its own pot for the duration of the 
experiment, making it possible to move the sample and reposition it accurately 

3) The samples also needed to be symmetrical about a centre of rotation. 

Equipment 
Detector 

The experiments were performed in the City University Radiation Laboratory using a Pantak 
HF160 industrial x-ray tube. An HPGe detector was used. This was because a good energy 
resolution was required for this experiment, to enable the resolution of the Compton and 
coherent peaks. A peak measured with the Ortec GLP-25300 HPGe detector, which was the 
detector used throughout all experiments, using an Am-241 source is shown in Figure 1 1 . The 
energy resolution is calculated as the FWHM of the peak, as illustrated in figure 11. For this 
detector the energy resolution is 0.435 keV at 59.54 keV (0.73%). The resolution for a Nal(Li) 
detector at the same energy is about 6-7%. The reason this peak was used to find the energy 
resolution is because the experiments that were carried out were done using the 57.97keV 
Kq2 peak from tungsten. This is very close in energy to the Am-241 peak at 59.54 keV and so 
the resolution will be approximately the same. 

■ 

Electronics 

A diagram of the electronics chain is shown in Figure 12. 

The detector was connected via a pre-amp and an amplifier to two single channel analysers, 
one to record the Compton peak and one to record a background region. Communication with 
a PC was enabled via an Ethernet card. 

Windowing 

An observed scatter spectrum is shown in Figure 1 3. The two coherent peaks of the Ka, and 
Kq2 W lines can be identified. The two, smaller Compton scatter peaks can be seen. 

The Ka 2 Compton peak was measured for this experiment. This is because the peak, 
although it has a stronger signal, is significantly overlapped by the two coherent peaks. 
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Transmission measurements 
Method 

The transmission measurements are a measure of the reduction in intensity of the 
unscattered peak and are a measure of the loss of counts due to tissue attenuation. For these 
measurements the detector was placed at zero degrees (See Figure 14). 

System calibration 

As the composition of the tissues being measured is unknown, the electron density 
measurement system needed to be calibrated. This was done by measuring some 
substances with a known or calculatable electron density. 5 substances were chosen in order 
to produce a comprehensive calibration graph. The solutions chosen were water, iso- 
propanol, and solutions of potassium hydrogen phosphate K 2 HP0 4 . Water and propanol were 
chosen because they are readily available; easy to handle and have a known electron density 
that is close to that of tissue. K 2 HP0 4 was chosen because it contains elements similar to 
those found in cellular fluids and so is a good model for human tissue composition. The 
concentration of the phosphate solutions could easily be varied to provide solutions with 
differing electron densities. In order to have values close to that of tissue, solutions of 2%, 5% 

■ 

and 10% were used. 

In order to verify the scatter data for the calibration solutions the linear differential scattering 
coefficient can be calculated theoretically as the composition of these solutions is known. 

The Klein-Nishina cross-section is dependent on incident photon energy and scattering angle. 
The Klein-Nishina differential scattering cross section is calculated to be 7.177 x 10* 29 
cm 2 /electron for 57.97 keV photons at a 28.2° scattering angle. Using this value and tabulated 
values for S(x) taken from Hubbell et al (1975) a value for Pcompton for each calibration solution 
was calculated. 

A graph showing the experimental values against the theoretical values is shown in Figure 15. 
The corrected scatter counts are the counts measured in the scatter peak corrected for 
attenuation and are given by 
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where Sco* is the counts recorded in the scatter peak corrected for attenuation. Smeas is the 
number of counts in the scatter peak, B s is the background counts in the scatter peak, l^as is 
the number of counts in the transmitted peak, B T is the number of background counts in the 
transmitted peak, l 0 is the unattenuated count intenstiy and B 0 is the background area for 
these counts. 

The graph of Figure 1 5 can be used to convert the corrected counts measured into differential 
scatter coefficients for Compton scatter,|j Sf where 



In the above equation (4.14) Scon- is the corrected scatter counts as described in equation 
4.13, N is the systematic experimental error and k is a constant that is found using the 
calibration graph. 

■ 

The trend line of the graph does not pass through zero but crosses the y-axis. This suggests 
that there is a systematic experimental error causing fewer counts to be recorded than 
expected. This is most likely due to a small amount of copper placed in the beam during the 
transmission measurements to protect the detector from the high photon flux. The geometry 
of the set-up was also changed between the scatter and transmission measurements. The 
detecor was moved further away. Due to the inverse square law this would mean that fewer 
counts would be recorded than expected. These two factors were not corrected for, as they 
can now be taken into account in this calculation. 

As the composition of the calibration solutions are known the electron v.3nsities of these 
solutions can be calculated using the following formula 



where p is the physical density of the material and ZJA is the ratio of atomic number to atomic 
weight for each element of mass fraction co. Z/A vales are tabulated and were taken from Attix 
(1996). The graph in Figure 16 shows the theoretical electron densities plotted against the 
measured scattering coefficients 

For this graph it can be seen that the two quantities correlate almost perfectly with a gradient 
equal to the Klein-Nishina cross section. This is what is expected as for high values of x the 
incoherent scattering factors become equal to Z. This agreement confirms the theoretical 
validity of the experiment. 




+ N 



44.14) 
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Results 

Figure 17 shows the results that were obtained from the scatter peak measurements. 

On the chart in Figure 17 the median of each tissue type is shown (thick middle line). The 
interquartile range is contained within the box and the whiskers show the total range. 

Analysis 

Calculation of electron density values 

The graph of Figure 1 5 gives a calibration equation to convert the number of counts in the 
scatter peak into the differential linear scatter coefficient p s . The equation given by the graph 
is 

ju s =1.737xl(r 7 jt + 7.919xl(T 3 44.16) 
where x is the corrected counts in the Compton peak. 

These experimental scatter coefficients are then converted into electron densities using the 
calibration solution values. This conversion is shown by the trend line in the results graph 
(Figure 1 1 ). 

The results ars. shown in the graph in Figure 1 1 . On this graph the values of electron density 
for standard tissue compositions given in ICRU report 44 (ICRU, 1989) are also displayed. In 
this report three separate values are given for different tissue compositions. The elemental 
compositions of these six tissues have been given in the table below. It is worth noting that 
the values quoted in this report are for healthy tissues only, as there is no published data for 
malignant tissue growths. 
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77?e elemental compositions (percentage by mass) of adult tissues 

(ICRU Report 44, 1989) 

It is usually assumed that malignant tissue has approximately the same structure as healthy 
glandular tissue. This is because tumours are usually within fibrous tissue rather than growing 
in purely fatty (adipose) tissue. 

The final results obtained are displayed in the table below. 



Tissue 


Electron density (e/cm 3 ) 


Benign 
Malignant 
Adipose 


(3.362 ± 0.141) x10 23 

(3.510±0.147)x10 23 

(S.S^iO.ISQJxIO 23 



Experimental values obtained for tissue electron densities 

This Difference in measured electron density between tissue types can be used in a model, 
such as the one described in our co-pending UK patent application GB0328870.1 , to 
determine the tissue type of samples for which the type is unknown. It therefore represents a 
potentially usefully diagnostic tool. As Compton scatter measurements can also be made in 
vivo, this approach also potentially lends itself to use as an in vivo, as well as in vitro, 
diagnostic approach. 

Although the three tissue types used to exemplify the invention here are •benign 1 , 'malignant' 
and adipose', the approach can be applied to the determination of other tissue characteristics 
or other tissue analysis applications. 
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EXPERIMENT B 
Material and Methods 

Samples 

A sample set of four different tissue types were examined comprising of 5 fibroadenoma 
(benign), 8 invasive ductal carcinomas (malignant), 4 fibrocystic change (non-malignant 
abnormal) and 5 pure adipose (normal) samples. Each sample was examined at two points. 
The samples were placed into polythene sample vials of 6mm inside diameter and 1mm wall 
thickness. Although the walls of the vial were relatively thick and would cause some 
attenuation of the scattered beam, these containers were chosen because they offered a 
number of important advantages. 

The sides were completely rigid so the samples could be placed into a vial and lightly 
compressed with a stopper without it distorting. This stopper is to remove any air gaps and it 
also minimises tissue movement throughout the experiment. The containers were cheap so 
each sample could have its own holder for the duration of the experiment, making it possible 
to move the sample and reposition it accurately. The samples also needed to be symmetrical 
about a centre of rotation. 

■ 

Method 

The Ka characteristic lines produced by a tungsten target x-ray tube were utilized as a 
monoenergetic source to ensure that the Compton scatter peak was detectable. Using this 
method the Compton and coherent scattered peaks from a recorded spectrum can be easily 
resolved and windowed and the bremsstrahlung background subtracted. The desired 
outcome of the experiment was to be able to resolve the Compton and coherent scattered 
peaks, whilst keeping them as close in energy as possible. The detector characteristics 
dictated that the minimum resolvable energy is about 1 keV. 

The experimental set-up is shown in figure 9. The x-ray beam was collimated to 0.5mm 
diameter, both before and after the sample. This was the smallest beam size viable whilst 
maintaining a reasonable flux. The !<„ line from the tungsten target of the x-ray source 
(Ete2=57.97 keV) was used. At this energy a scattering angle of 30° gave a peak separation of 
1 keV between the Compton and coherent scatter peaks. The scattering volume comprises of 
the tissue contained within the intersecting area of the incident and scattered beam. For this 
beam collimation and scattering angle the entire scattering volume was contained within the 
sample, with no air or polythene from the vial included. Each sample was measured for a total 
time of four hours, with the sample being rotated throughout the measurement in order to 
reduce any errors due to the inhomogeneity of the tissues. 
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Equipment 



The experiments were performed using a Pantak HF160 industrial x-ray tube. An HPGe 
detector was used in order to produce the energy resolution required to resolve the Compton 
and coherent peaks. The energy resolution was measured to be 0.435 keV at 59.54 keV 
(0.73%). The detector was connected via a pre-amp and an amplifier to two single channel 
analysers, one to record the Compton peak and one to record a background region. An 
observed scatter spectrum of a malignant tissue is shown in figure 20. 

The two coherent peaks of the Kq, and Ko 2 W lines can be identified and the two smaller 
Compton scatter peaks can be seen. The Kcc Compton peak was windowed over an area 
where there was no superposition of the coherent peak. This windowed area, which was 
used for the scatter measurements, is also shown in figure 20. The transmission 
measurements for each sample were made by placing the detector at zero degrees and 
recording the photon intensity with and without a sample in position in the beam. 



System calibration 



As the composition of the tissues being measured is unknown, the electron density 
measurement system needed to be calibrated. This was carried out by measuring substances 
with a known electron density or one that could be calculated. Five substances were chosen 
in order to produce a calibration curve. 

The solutions chosen were water, iso-propanol, and solutions of potassium hydrogen 
phosphate K 2 HP0 4 . Water and propanol were chosen because they are readily available, 
easy to handle and have a known electron density that is close to that of biological materials. 
The concentration of the phosphate solutions could be varied to provide solutions with 
differing electron densities. In order to have values close to that of tissue, solutions of 2%, 5% 
and 10% were used. 

In order to verify the scatter data for the calibration solutions the linear differential scattering 
coefficient can be calculated theoretically as the composition of these solutions is known. The 
linear scattering coefficient is a measure of the probability that a photon of incident energy E 
will be scattered through an angle 9 and is given by equation (7): 



/W~=P^^f (7) 
WHere i = Z^, (8) 
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where M is the molecular mass of the material, p is the mass density; and N A is Avogadro's 
constant. S(x) is the incoherent scattering factor and the differential scattering cross section is 
denoted with KN for the Klein-Nishina cross section. The Klein-Nishina differential scattering 
cross section for the Compton effect is given by 



da 



KN _ ' 



1 



13 



dQ Ll + a(l-cos#) 



1 + COS0 



, or (l- cos #) 2 

{l + cos 2 0)[l + a(l-cos0)] 



(9) 



where E is the incident photon energy and 9 is the photon scattering angle, a is the ratio of 
the incident photon energy to the electron rest mass energy given by 



a = 



m 0 c 



(10) 



and r 0 is the classical electron radius. 

The Klein-Nishina differential scattering cross section is dependent on photon energy and the 
angle of scatter. It was calculated to be 7.177 x 10' 26 cm 2 /electron for 57.97 keV photons at a 
30° scattering angle. Using this value and tabulated values for S(x) taken from Hubbell et a/. 
(1975) a value for Mcompton for each calibration solution was calculated. A graph showing the 
experimental scatter measurements against the scatter coefficient values calculated from 
equation (7) is shown in figure 21. 

The corrected scatter counts are the counts measured in the scatter peak corrected for 
attenuation and are given by 



^corr 



[ S mcas ~ B s ] 



I mens 



(11) 



where is the counts recorded in the scatter peak corrected for attenuation. is the 
number of counts in the scatter peak. B s is the background counts in the scatter peak, l meas is 
the number of counts in the transmitted peak, Br is the number of background counts in the 
transmitted peak. I 0 is the unattenuated count intensity and B 0 is the background area for 
these counts. Figure 21 can be used to convert the corrected counts measured into 
differential scatter coefficients for Compton scatter,p s> where 



M,=k[S„„] + N 



(12) 
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In equation (12) is the corrected scatter counts as described in equation (1 1), N is the 
systematic experimental error and k is a constant that is found using the calibration curve. 
As the composition of the calibration solutions are known the electron densities of these 
solutions can be calculated using the following formula 



p.-p"jLrfr* (13) 



where p is the physical density of the material and Z/A is the ratio of atomic number to atomic 
weight for each element of mass fraction w. Z/A values are tabulated and were taken from 
(Attix 1986). Figure 22 shows the theoretical electron densities calculated from equation (13) 
plotted against the measured scattering coefficients given by equation (12). It can be seen 
that the two quantities correlate well with a gradient equal to the Klein-Nishina cross section 
which is expected, as for high values of x the incoherent scattering factors become equal to Z. 

Results 

Figure 23 shows the results of the electron density measurements that were obtained from 
two points on each sample. The median of each tissue type is shown (thick middle line). The 
interquartile range is contained within the box and the whiskers show the total range. 
The graph of figure 21 gives a calibration equation to convert the number of counts in the 
corrected scatter peak into the differential linear scatter coefficient u s . The equation given by 
the graph is 



M, =1.737xl0" 7 x + 7.919xl0" 3 (14) 
where x is the corrected counts in the Compton peak. 

These experimental scatter coefficients are then converted into electron densities using the 
Klein-Nishina cross section. This conversion is shown by the trend line in figure 22. The 
average results are shown in figure 24. 

In figure 24 the values of electron density for standard tissue compositions given in ICRU 
report 44 (ICRU 1989) are also displayed. In this report three separate values are given for 
different tissue compositions. The elemental compositions of these six tissues have been 
given in table 4. It is worth noting that the values quoted in this report are for healthy tissues 
only, as there is no published data for malignant tissue growths. 
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Table 4. The elemental compositions (percentage by mass) of adult tissues (ICRU Report 
44, 1989) 



It is usually assumed that malignant tissue has approximately the same structure as healthy 
glandular tissue. This is because tumours are usually within fibrous tissue rather than growing 
in purely fatty (adipose) tissue. 
The final results obtained are displayed in table 5. 



Tissue 


Electron density (e/cm 3 ) 


Benign 


(3.330 ± 0.140) x10 23 


Malignant 


(3.490±0.147)x10 23 


Adipose 


(32Q1±0A3Q)xW 2Z 


Fibrocystic 




change 


(3.752±0.158)x10 23 



Table 5. Experimental values obtained for tissue electron densities 

Each individual measurement is subject to statistical variation. The error o is given 




where x is the mean number of counts if the reading is repeated N times. For 



the scatter readings each measurement was measured for a sufficient time (4 hours) to 
ensure that the error on the counts was sufficiently low (<0.5%). Due to time constraints the 
readings were not repeated. 

The largest error is associated with the subtraction of the background counts. The overall 
error on the background count calculation is 4.2%. This is shown by the error bars in figure 
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21 . Other sources of error are the effect of multiple scatter, the error in positioning and the 
error in repositioning the sample for the transmission measurements. There is also a widening 
of the Compton scatter peak caused by the acceptance angle of the pre-detector collimator. 
None of these other errors have been considered as they are difficult to quantify and are 
small compared to the background subtraction error outlined above. 

Discussion and Conclusions 

The results show that there is a detectable difference between the electron density of adipose 
and malignant tissue, to a value of 6.4%. This difference is consistent with the values found 
by using the adipose and glandular tissue values from ICRU report 44. The average value for 
glandular tissue (ICRU44 glandular#2) is 6.2% higher than the average adipose value 
(ICRU44 adipose #2). 

There has been no composition values published for benign (fibroadenoma) or fibrocystic 
tissues. However the measurements made within this experiment B found a difference In the 
electron density of benign and malignant tissues to the value of 5.6% and a difference 
between fibrocystic change and malignant tissue to be 2.3%. 

It is difficult to verify these results using the literature as there is no published data from any 
previous studies using this tissue type. However the high degree of correlation for the 
calibration solutions (figure 21) shows that the system has a reliable accuracy. This inspires 
confidence in the findings that there is a measurable difference between the benign and 
malignant tissues. 

There is a great deal of evidence to suggest that the metabolism and physiology of tumour 
cells differ greatly to that of normal and indeed benign cells. 

Within a benign tumour growth there is often an increase in cell proliferation but the cells 
themselves are relatively normal. However, in a malignant lesion the structure and 
metabolism of the tumour cells and host tissue have a different biochemical structure (Gould 
1997). This implies that the increase in the electron density of benign tissues compare to 
normal may potentially be due to an increase in cell concentration rather than to changes in 
composition, as seen in malignant tissues. This is consistent with the finding that benign 
tissues ex-vivo have an electron density which is only slightly higher than normal tissues and 
malignant tissues display a much larger difference. 

Dr Otto Warburg first observed in 1930 that cancer cells have a fundamentally different 
energy metabolism than normal cells (Warberg 1930). Since then research has shown that 
tumour cells undergo anaerobic glycolysis, the process where glucose is converted to lactic 
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acid through the process of fermentation. This process is extremely inefficient compared to 
normal cell aerobic respiration. 

Glucose consumption rate has been shown to be proportional to histological grade (Vaupel et 
al. 1989) and high grade tumours can absorb about 40 times more glucose in order to supply 
their high energy demands for increased growth. This process is what makes positron 
emission tomography (PET) imaging so effective at imaging tumours using 18F-FDG, an 
analogue of glucose. It enables PET to distinguish between benign and malignant neoplasms 
with a high degree of accuracy, as benign tissues do not exhibit increased glucose 
consumption (Brock et at. 1 997). 

Anaerobic glycolysis causes a build up of lactic acid to occur within the tissue. The lactic acid 
(CH3-CH(OH)-CO(OH)) which builds up within the tumour has a high electron density 
compared to the host tissue of 8.2 x 10 23 electrons/cm 3 and so could be responsible for the 
increase in electron density that is measured. There is also an increase of ketones and 
glutamine (Vaupel et al. 1 989) which may also increase the overall electron density of tumour 
tissues. Although no direct measurements have been made of the composition of benign and 
malignant tissues, the above suggests that there are significant differences in composition It 
is difficult to estimate the precise nature of the composition changes, given that there are a 
number of processes occurring in the tissue during tumourgenesis. 

The final tissue type that was examined was fibrocystic change. Although this term 
encompasses a range of histological changes, the majority are characterised by tissue 
fibrosis. This is a scarring process whereby the stromal (connective tissue) component of the 
tissue is increased and collagen accumulates. Although increased mature collagen may be 
seen in a few other benign disease processes in the breast, the most pronounced increase 
probably occurs during fibrocystic change. This may account for the finding that this tissue 
classification had a higher electron density than any other W a of tissue, even malignancy 
When examining the tissues exhibiting fibrocystic change it is likely that any fluid filled 
pockets (cysts) will become dispersed during tissue preparation leaving only the dense fibrotic 
tissue under examination. 

Not only is the present invention useful for assessing increased numbers of fibroadenomas 
.nvasive ductal carcinomas and FCC tissues, the present invention may also be adapted to' 
assess healthy fibrous tissue and further disease processes. 
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